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Rolling and sliding of a nanorod between two planes: Tribological regimes and control 

of friction 

Mykhaylo Evstignee\G and Peter Reimann 
Universitat Bielefeld, Fakultat fur Physik, 33615 Bielefeld, Germany 

The motion of a cylindrical crystalline nanoparticle sandwiched between two crystalline planes, 
one stationary and the other pulled at a constant velocity and pressed down by a normal load, is 
considered theoretically using a planar model. The results of our model calculations show that, 
depending on load and velocity, the nanoparticle can be either rolling or sliding. At sufficiently high 
normal loads, several sliding states characterized by different friction forces can coexist, correspond- 
ing to different orientations of the nanoparticle, and allowing one to have low or high friction at the 
same pulling velocity and normal load. 

PACS numbers: 46.55.+d, 62.20.Qp, 68.35. Af, 07.79.Sp 
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I. INTRODUCTION 

Because of its great importance for nanotechnological 
applications, control of friction at the nanoscale is a hot 
topic of current research. While conventional lubricants 
cannot be applied for this purpose, research efforts have 
been guided by the vision of Richard Feynman that the 
nanobearings can "run dry"—. It has been shown that 
decreasing the normal load 2 or pulling velocity^, as well 
as normal load actuation^ can lead to a dramatic fric- 
tion reduction. A mechanism of friction control par- 
ticularly pertinent to our present work is the so-called 
structural lubricity, or superlubricitySr— , arising due to 
the structural incommensurability of the two contact- 
ing surfaces. More precisely, each atom of the sliding 
surface feels the force generated by the periodically ar- 
ranged atoms of the substrate, so that the total friction 
force is the sum of the forces felt by each atom of the 
slider. If the slider and the substrate are incommensu- 
rate, these forces add up randomly, resulting in nearly 
frictionless sliding. The phenomenon of superlubricity 
has been observed experimentally^ 1 ^— . On the other 
hand, it has been demonstrated by means of molecular 
dynamics simulations 14 , as well as stochastic modelling 
and experiment— that a flat nanoobject (e.g. a graphite 
flake) in contact with the surface quickly reorients it- 
self into the "commensurate" state of high friction, even 
though the initial orientation may be the superlubric one. 

Finally, one can influence the friction forces by using 
the rolling motion of round nanoparticles in between the 
surfaces 1 ^—. Nanoparticle rolling can be identified ex- 
perimentally by making a small indentation on its surface 
with a sharp tip of an atomic force microscope and subse- 
quent localization of this mark after the manipulation 19 . 
Alternatively, for a highly symmetric nanoparticle, such 
as a carbon nanotube, a rolling regime can be identified 
from the characteristic periodicity of the time-dependent 
friction force - in the rolling state, this periodicity is pro- 
portional to the circumference of the nanoparticl o 20 i 21 . 

In this paper, we consider the motion of a cylindrical 
nanoparticle sandwiched between two planes, see Fig. Q] 
Because of the unavoidable deviations from a perfect ro- 
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FIG. f: Schematic illustration of the system: a nanoparticle is 
sandwiched between the stationary lower plane and the upper 
plane pressed down by the normal load /jv and attached to 
a spring of stiffness k, whose other end is pulled with the 
velocity V. The friction force / is deduced from the elastic 
deformation of the spring attached to the upper plane, cf. 
Eq. ©. 



tational symmetry, the surface of the nanoparticle ac- 
tually consists of facets characterized by different lat- 
tice constants, and hence by different commensurabilities 
with both planes^ 1 -. Based on a simple model described 
in Section HH we show in Section Mil that the nanopar- 
ticle can be stabilized in several friction states, namely, 
rolling friction and different sliding regimes that can be 
realized depending on which facets of the nanoparticle 
are in contact with the upper and lower surface. These 
results are summarized in a state diagram showing the 
stability regions of different friction states depending on 
the two experimental control parameters - normal load 
and pulling velocity. The state diagram contains a region, 
where sliding friction regimes characterized by commen- 
surate contact (high friction forces) and incommensurate 
contact (low friction forces) of the nanoparticle coexist. 
We propose a friction switching scenario, allowing one 
to realize either low or high friction forces at the same 



pulling velocity and normal load. 



II. THE MODEL 
A. Equations of motion 

We consider a crystalline nanoparticle, approximating 
as close as possible a round shape of some preset radius, 
see Fig. [TJ Unlike the symmetric fullerene-like molecules 
with equivalent 'facets' studied in Refs. llftlla . in our sys- 
tem, the facets have different number of atoms and differ- 
ent commensurability with the planes, see Fig.Q] so that, 
depending on which facets are in contact with the planes, 
different friction regimes may be realized. While our 
model system from Fig. [T] is two-dimensional, the main 
results immediately carry over to some three-dimensional 
objects, such as nanorods. We henceforth consider a sin- 
gle particle, expecting analogous findings for the case of 
many particles, provided they are sufficiently dilute. 

The overall geometric configuration of the system is 
completely determined by the nanoparticle's center of 
mass r = xe x + ye y , its rotation angle 0, and the po- 
sition R = Xe x + Ye y of some reference atom from the 
top plane, which we assume to be horizontal at all times; 
here e x . y are the unit vectors in the x- and y-direction. If 
the nanoparticle and the planes could be viewed as rigid 
bodies, the equations of motion for the generalized coor- 
dinates (r, <j), R) would be derivable from the Lagrangian 

w , x, • i ^x rnr 2 ift MR 2 TT , ,. 
L(r,&R,r,^R) = — + -|- + — U B (r,4>) 



-U T{ R-r,t)- K -^£ + f N Y, 



(1) 



where the first two terms represent the translational and 
rotational kinetic energy of the nanoparticle with mass 
m and moment of inertia /; the third term describes the 
kinetic energy of the upper plane of mass M; the fourth 
and the fifth terms correspond to the interaction between 
the nanoparticle and the bottom and top planes, respec- 
tively; the last two terms correspond to the energy of 
elastic deformation of the spring of stiffness k whose other 
end is pulled at the constant velocity V, and the energy 
of the upper plane produced by the normal load. There- 
fore, the elastic force generated by the spring is equal in 
magnitude to the instantaneous friction force 



/ = -k(X - Vt) 



(2) 



We furthermore assume that the separation between the 
two planes - the nanoparticle's diameter - is sufficiently 
large, so that the interaction energy between them is neg- 
ligible. 

In reality, the atoms of the nanoparticle and the planes 
are not rigidly coupled to each other, and their motion af- 
fects the motion of the global coordinates r, 0, and R. If 
the time scales of the overall nanoparticle motion is much 



slower than the time scale of individual atoms, the inter- 
action between the global degrees of freedom, r, </>, R, and 
those of atoms composing the nanoparticle, both planes, 
and the spring can be approximately taken into account 
by means of the following three modifications of the equa- 
tions of motion generated by the Lagrangian (flj 22 : 

(i) renormalization of the forces and torques acting on 
the relevant coordinates r, (f>, and R; 

(ii) introduction of velocity-dependent dissipation 
forces describing the effect of energy loss from the global 
degrees of freedom r, <fi, and R into the atomistic de- 
grees of freedom of the nanoparticle, both planes, and 
the spring attached to the upper plane; 

(iii) introduction of noise corresponding to the inverse 
process of energy transfer from random atomic vibrations 
into the global degrees of freedom. 

The effect (i) is accounted for by reinterpreting the en- 
ergies Ub,t as free energies. The dissipative forces (ii) 
due to the internal degrees of freedom of the nanopar- 
ticle and the planes can be derived from the dissipation 
function^, 



$(r, (j), R, f , <£, R) = - ^2 Vqg' (r, </>, R) qq' 

9:9' 



(3) 



where the generalized coordinate indices q, q' run over 
the values (x, y, (p, A, Y), and the dissipation coefficients 
are symmetric, 



Tjqq' — Tjq'q 



(4) 



These forces are written as velocity derivatives of the 
dissipation function, /* 5S = —d&/dq = — ^q'Vqq'q'- 
Similarly, an additional dissipative force on the upper 
plane arises due to the internal degrees of freedom of the 
spring: ff^. 8 ing = — ?7s(R— Ve x ). This additional dissipa- 
tion channel is not present in the dissipation function ([3]) , 
because the position Vt of the pulled end of the spring 
is not a generalized coordinate in the Lagrangian (fl]); 
therefore, the dissipative force f^lng should be included 
"by hand". Finally, we account for the noise effect (iii) 
by adding suitably chosen Gaussian white noises to the 
right-hand side of the equations of motion: 

mr = -Vp7 B (r, 0) - I7 T (R - r, 0)] 

-■q rT T - 77 r 00 - ?7rRR + £r(£), 

fill 

I<f>-. 



d[U B (v, ( j ) ) + U T (R-r, ( j } )} 



MR = -W T (R - r, </>) - k(X - Vt)e x - f N e y 
-VRrr ~ 77R00 - 7/ RR R - r)s(R - Ve x ) 

+Ca(t) + &(*)• (5) 

Here, rj rr is a tensor with components r] xx ,rj xy ,r/ yx = 
TixyiVvyi with similar definitions for t^r and ?)rr. Also, 
r] r4> and ?7r,0 are vectors, e.g. r) r( p = r\ x ^& x + rj y 4,e y with a 
similar definition for rj R( p, and -q^ = rjj^, i]^ = 77^ 



are the respective transposed vectors. The Gaussian 
white noises £, q (t), q — x,y,<fi,X,Y have zero mean and 
obey the fluctuation-dissipation theorem of the second 
kind, (£g(*)£«2'0')> = 2Tr] qql (r,^,R)S(f -t). The noise 
£<?(£) due to the spring is uncorrelated with the noises 
£q(t) and its autocorrelation function is (£sa(t)£sp(t')) = 
2Trjs5 a [i5{t' — t), where a, ft refer to the coordinates x, y. 



B. Choice of the functional form for the potentials 
and dissipation coefficients 

Having written down the general equations of motion 
([5]) for the relevant coordinates, we should specify the 
functional forms of the potentials Ub,t, as well as the 
dissipation coefficients r] rr , r] r( p, etc. We consider the po- 
tentials first. For simplicity, we assume the upper and 
the lower planes to be equivalent, leading to 

U T (R-r,cf> + ir) = U B (r,(f>)=:U(r,cf>) , (6) 

and leaving us with the necessity to determine the func- 
tion U(r,(j>). This function should possess the symme- 
tries 

U(t, <j>) = U(r + ae x ,<f>) = U(r, (f, + 2n/N) , (7) 

where a is the lattice constant of the plane, and N is an 
integer related to the rotational symmetry of the parti- 
cle. For example, for a nanoparticle as sketched in Fig.[TJ 
we have N = 4. There are many functions with the 
property ([7]). To make a physically motivated choice, 
we view the nanoparticle as a collection of periodically 
arranged "pseudoatoms" , where each pseudoatom repre- 
sents a group of closely arranged real atoms of the par- 
ticle. In particular, if Fig. [1] depicts a cross-section of 
a nanorod, then each pseudoatom describes the cumula- 
tive effect of an atomic row along the axis perpendicular 
to the plane of the figure. The ith pseudoatom of the 
nanoparticle has the coordinate 



r;(r,( 



■diei((f>) , 



(8) 



where di = |r» — r| is the rigidly fixed distance of the 
ith pseudoatom from the center of mass, and ej(0) = 
e x cos(9i + </>)+ e y sin(#j + <j)) is the unit vector pointing 
from the center of mass to the ith pseudoatom. The con- 
stant angles 9i are completely determined by the crystal 
structure of the nanoparticle. 

The potential of interaction with the bottom plane is 
the sum of the respective interaction energies of all pseu- 
doatoms 



J7(r,0=53u(ri(r,$) 



(9) 



The functions u(r) = u{v + ae x ) can be expanded into 
Fourier series in x with the expansion coefficients depend- 
ing on the 77-component of r. Neglecting the second and 
higher harmonics, we employ the functional form 

u(r) = u Q (y) + ui(y) cos . (10) 



To account for the possibility of adhesion, the zero-order 
term is taken to be the Lennard- Jones potential, 



MV) = e {(a/y) 12 - 2{a/yf) 



(11) 



where e is the adhesion energy and a the equilibrium 
separation from the surface. 

The function U\{y) in Eq. (fTU|) has the physical mean- 
ing of the corrugation amplitude of the potential (TTOl) in 
the x-direction. We assume it to increase exponentially 
upon approaching the surface: 



ui(y) = AUe- {v - a)/x , 



(12) 



where AC/ is the corrugation at the equilibrium separa- 
tion cr, and A is the characteristic decay length. 

With respect to the dissipation coefficients of the 
nanoparticle, we should distinguish between the contri- 
butions due to the bottom and the top planes. We as- 
sume that in the course of its translational and rotational 
motion, the total dissipative force on the nanoparticle is 
a sum of the respective contributions from all its pseu- 
doatoms. That is, the motion of the ith pseudoatom with 
the velocity r^ relative to the bottom plane results in the 
dissipative force —rj{vi)Yi on that pseudoatom. The ex- 
pression for the dissipative force contribution from the 
top plane is similar, but with the dissipation coefficient 
77 (R — Ti) and the relative velocity r.^ — R. 

For simplicity, we assume isotropy of the damping co- 
efficients r)(ri), which are treated as scalar functions of 
the position rj. We choose the following functional form: 



17(1-0 = rjoe-to-'M 



(13) 



770 being the damping coefficient at the minimum of the 
Lennard- Jones potential (|TT|) . and £ the decay length. 
Explicitly, the velocity of ith pseudoatom is expressed in 
terms of the generalized velocities r and <fr as 



r*(r,< 



r + dit, 



(14) 



where the tangential vector is 



u(d>) 



dei{4>) 



= -e x sm(6i + 4>)+e y cos(9 t + (j>) . (15) 



The dissipative force due to the bottom plane is 

i i i 

. (16) 

The dissipative force due to the top plane has a similar 
form, but with the velocity r^ replaced with r^ — R and 
77(1*0 with 7y(R — i*0, that is 

fL s = -^-7(R- r *)(r-R)-I]^R-r0rfiW- (17) 

i i 

Summing all the contributions, we obtain: 
?7rr = J2[v(*i) + ri(K - i*0]I , 



i 

r)rR = -J2n{R-ri)I, (18) 

i 

where I is a unit 2x2 tensor. 

In view of the symmetry Q of the dissipation co- 
efficients, the second of these equations (|18[1 uniquely 
fixes rj^r = rjj,. To determine the remaining coefficients 
7/0^ and ?70R, we should consider the torque produced 
by the dissipative force as the nanoparticle rotates be- 
tween both planes. The magnitude of this torque is 



Kdis 



Y^idi r](r t )r l +7](R-r l )(r i 



R) 



substitution of the expression (TlH) we find 



•tj. Upon 



V<t>R = 






Tl(R-Ti)ditl 



[^(rO+^R-r,-)]^ 



(19) 



For the upper plane, we have from the symmetry (U) 
of the dissipation coefficients ?^R. r = r; r R and ?7r,0 = ??T R . 
The remaining tensor 77RR, describes the effect of energy 
dissipation of the upper plane into the internal degrees of 
freedom of the stationary nanoparticle and has the form: 



»)RR = Xl 7? ( R_r ^ 1 



(20) 



C. Overdamped zero-temperature limit 

It is difficult to estimate the dissipation coefficients 
from first principles, because the basic building block of 
our model - a pseudoatom - is a complex object con- 
sisting of many real atoms. It is not unreasonable to 
assume though that the dissipation coefficient of such a 
pseudoatom (and of the nanoparticle itself) can be many 
orders of magnitude higher than that of a true atom 
on a surface. Therefore, in our numerical calculations, 
we assume that the dissipation effects are much stronger 
than the inertia effects, allowing us to consider the over- 
damped limit by formally setting the nanoparticle's mass 
and moment of inertia to zero: m = 0, J = 0. Likewise, 
we assume that the spring attached to the upper plane 
is overdamped, allowing us to set M = 0. Finally, since 
the potential energies from Eq. (TTJ) represent an effect 
of many atoms, noise effects can be assumed extremely 
small in comparison to the interaction forces and the nor- 
mal load. Therefore, we neglect thermal noise by setting 
T to zero in the equations of motion (O, yielding 

77 rr r + r) T $<fi + 77 rR R = 

-V[[/(r,0)-[/(R-r> + 7r)] , 

%> r r + r/^ip + ^rR = 

d[U(r, 4>) + U(H - r, <f> + ir)] 



?7R r r + VR^ + ?7rrR = -VC/(R - r, cj) + n) 

-r) S (R - Ve x ) - k(X - Vt)e x - f N e v . (21) 

These equations can be simplified even further if we 
consider a nanoparticle, which is symmetric with respect 
to rotations by n, as in Fig. [TJ Due to this symmetry, and 
due to the equivalence of the upper and lower planes, we 
can state that there is a solution of Eqs. (|2"Tj) . for which 



r = R/2 



(22) 



up to an addition of an integer multiple of the lattice 
constant a in the x-direction. This is verified by inspec- 
tion of the first of the equations of motion (|21[) . where 
substitution of the relation (|2"2")l renders the force in the 
right-hand side vanish. Considering the left-hand side, 
let us have a closer look at the damping coefficients 
from Eq. (|18[) . Because our nanoparticle is symmetric 
with respect to rotations by w, for each pseudoatom at 
Yi = R/2 + diei there is a symmetric partner pseudoatom 
at Tk = R/2 — diei. Then, comparison of the first and the 
third equations (|18[) yields r) rr = — 2?7 r R. Furthermore, 
since the tangential vectors (fT5|) of the ith and the fcth 
pseudoatoms are opposite to each other, tj, = — tj, the 
sum in the second equation (fTHJ) vanishes, r] r $ = 0, auto- 
matically implying that r)^, r — 0. Then, the first equation 
(I2T1) reduces to ry rR (R - 2r) = 0, implying Eq. (|22J}. By 
numerically integrating the full set of equations (|2T|) , we 
have verified that the relation (|22p is stable: for all ini- 
tial conditions tried, the system eventually entered the 
regime with r = R/2. Therefore, computational effort 
can be reduced by roughly a factor of two by replacing 
five equations (|2ip with three equations of motion for 
(6, X, and Y: 



%,rR 



dU(R/2, 



[VS + %R + W2]R + VR<f><t> = - W(R/2, 4>) 
+ [ m V-K(X - Vt)]e x - f N e y . (23) 

III. RESULTS AND DISCUSSION 
A. Parameters and units 

In all our numerical results below, we have chosen the 
lattice constant a as the unit of length, the adhesion en- 
ergy £ as the unit of energy, and the ratio 7] a 2 je as the 
unit of time. This choice fixes the unit of force to e/a, 
the unit of velocity to e/(ar]o), and the unit of spring 
constant to e/a 2 . The value of the dissipation coefficient 
at the minimum of the potential (|11[) in these units is 
r/o = 1, and, obviously, a = 1, e = 1 in these units. 

Our nanoparticle is constructed from an arrangement 
of pseudoatoms in a square lattice with a lattice con- 
stant 6, which, in general, is not equal to the lattice 
constant a of the two planes. From this lattice we se- 
lect those pseudoatoms whose distance from the cen- 



ter of mass is smaller than the preset radius. In nu- 
merical simulations of Eq. (|23p . we focused on an ap- 
proximately round crystalline nanoparticle with a radius 
of 5b. The "commensurate" facets of the nanoparticle 
correspond to the rotational angle given by an integer 
multiple of n/2, 4> C omm — nir/2, whereas the "incom- 
mensurate" facets correspond to its half-integer multiple, 

4>incomm = (n + l/2)n/2. 

In our simulations, we have tried different values of b, 
and obtained qualitatively the same behaviour as for the 
b = a = 1 case reported below. As for other parame- 
ter values, we have taken the equilibrium distance of the 
Lennard-Jones potential (jlip to be equal the lattice con- 
stant, a — 1. The corrugation depth of the potential (|T2|) 
was taken to be AU = 3/4, and its decay length, as well 
as the decay length of the damping coefficient (TT4")) were 
set to A = £ = 1/5. Finally, the spring constant was set 
to k = 1, and the spring damping coefficient was r)$ = 10. 
Other values of these parameters of a comparable order 
of magnitude produced qualitatively similar results. 
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FIG. 2: State diagram of the system from Fig. Q] showing the 
stability regions of the low friction sliding (LFS), high friction 
sliding (HFS), and rolling (R) states. 



B. State diagram 



tion, which can be either the low-friction sliding (LFS) 
or the high-friction sliding (HFS) state. In both these 
states, the orientation angle <p of the nanoparticle per- 
forms small rocking motion around the value, which is 
either 7r/4 (LFS state) or (HFS state) , up to an integer 
multiple of tt/2. The LFS state corresponds to the in- 
commensurate contact between the nanoparticle and the 
planes, while the high-friction sliding (HFS) state corre- 
sponds to the commensurate contact. The LFS state is 
stable above the solid line in Fig. [3J the HFS state above 
the dashed line. Above both separation lines, both LFS 
and HFS states are stable, and the actual state of mo- 
tion of the nanoparticle depends on its initial prepara- 
tion. Interestingly, to the left of the intersection point 
at ca. V = 6, there is a region at lower normal loads 
where the LFS state is the only stable state of motion; 
similarly, to the right of this point, there is a region at 
higher normal loads where the particle can exist only in 
the HFS state. 

We now try to understand the state diagram from 
Fig. [5] In order to destroy the contact, work must be 
performed against the adhesion forces and the normal 
load. This wok is done by the viscous drag, which scales 
linearly with the velocity V, and by the force generated 
by the moving corrugated potential. At low pulling veloc- 
ities, it is the latter force that is responsible for breaking 
the contact. If the nanoparticle is in the LFS-state, the 
effective corrugation of the potential is lower than that 
in the HFS-state, because the contact is incommensurate 
and involves fewer atoms. Consequently, a smaller nor- 
mal load is required to stabilize the LFS state, explaining 
its stability island in the low- velocity region. 

At faster pulling, on the other hand, the effect of the 
potential corrugation becomes less important. This is so 
because the upper plane moves relatively fast with re- 
spect to the particle, so that the particle cannot follow 
the fast temporal variations of its potential and feels, 
instead of the true corrugation depth u\ , a smaller time- 
averaged corrugation. Therefore, it is the viscous drag 
that is responsible for turning the nanoparticle at fast 
pulling. In order to break the contact, one has to over- 
come the adhesion between the nanoparticle and both 
planes. Since adhesion forces are larger in the HFS-state, 
smaller normal load is required to stabilize this state at 
fast pulling. 



We have found that depending on normal load /at and 
pulling velocity V, the nanoparticle can be rotating, or 
it can be stabilized in a sliding state of either low or high 
friction. The results of our numerical simulations are 
summarized in the state diagram of the system, Fig. [5J 
showing which friction regimes are stable for given values 
of /jv and V. When the normal load is sufficiently low, 
the nanoparticle can only exist in the rolling (R) state 
of motion, where the magnitude of the orientation an- 
gle steadily grows in time, as in Fig. 0Jb). On the other 
hand, high normal load stabilizes the sliding state of mo- 



C. Friction force 

Fig. [3] exemplifies the typical evolution of the friction 
force ([2]) in the three states. The curves (a) and (b) re- 
semble the typical evolution of the friction force in an 
atomic friction experiment in the stick-slip regim o 24 ' 25 . 
During the stick phases, the nanoparticle and the upper 
plane are almost stationary, while the elastic deforma- 
tion of the spring attached to the top plane constantly 
increases due to pulling. When elastic energy becomes 
sufficient to initiate the slip, the nanoparticle gets dis- 
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FIG. 3: Temporal evolution of the friction force ((2| in (a) 
the LFS-state, (b) the HFS-state, and (c) the R-state. All 
curves are obtained for the same pulling velocity V = 0.1, 
but different normal loads: /jv = 100 for the curves (a) and 
(b), and /jv = 1 for the curve (c). The difference between the 
LFS- and HFS-curves (a) and (b) is in the orientation angle 
of the nanoparticle: namely, for the curve (a), the angle <j> is 
close to it/ 4, and for the curve (b), it is close to zero. 



placed in the ^-direction by one lattice constant, and the 
upper plane by two lattice constants, see Eq. (|2"2"|) . result- 
ing in a sudden relaxation of the spring. Consequently, 
the periodicity of the stick-slip curves in the LFS regime 
(a) and the HFS regime (b) equals two lattice constants. 
We note that the modulation amplitude of the stick-slip 
curves (a) and (b) is about the same, while the mean fric- 
tion forces developed in both regimes are very different, 
in spite of the fact the pulling velocity and the normal 
loads are identical for the curves (a) and (b). 

The curve (c) depicting the evolution of the friction 
force in the rolling state is interesting in three respects: 
its periodicity is notably larger than in the LFS- and 
HFS-cases (a) and (b); its modulation amplitude is also 
surprisingly high; the slope of the curve in a stick phase 
is much smaller than in the cases (a) and (b). These 
peculiarities can be explained as follows. 

The periodicity of the "stick-roll" curve (c) is equal 
to the distance travelled by the upper plane during the 
"slip"—. This distance is quite large, because the "slip" 



event, seen as the sudden drop of the friction force, is 
associated not with the transition of the nanoparticle by 
one lattice constant, but with its rotation by an angle 
7r/2 from one commensurate contact state to the next. 
During this rotation, the particle's center of mass travels 
a distance of about irR/2. The actual distance is slightly 
smaller than this value, because the nanoparticle is not 
perfectly round. In view of Eq. (|22p . the periodicity of 
the curve (c) twice that value, that is, slightly smaller 
nR. For R = 5, this gives the periodicity that should be 
a bit smaller than 16 lattice constants. The periodicity 
of the curve (c) is indeed 14 lattice constants. 

The second interesting feature is a roughly ten-fold 
larger amplitude of force variations in the rolling state 
compared to the LFS- and HFS-curves in Fig. [31 in spite 
of the fact that the normal load is two orders of mag- 
nitude smaller. As it turns out, the reason is precisely 
the much smaller normal load. In each stick phase, the 
elastic energy of the spring constantly builds up and is 
suddenly released to induce a rotation of the nanoparti- 
cle in the rolling state from Fig. [3] (c). At the same time, 
the nanoparticle in the stick phase can slightly turn and 
lift the top plane up against the small normal load. Since 
part of the torque applied to the particle by the spring is 
used to lift the upper plane, a much larger elastic force 
is necessary to induce the particle's rotation, implying a 
large modulation amplitude of the curve in Fig.^c). For 
the curves from Fig. [3] (a) and (b), on the other hand, 
the normal load is too high to allow for any significant 
lifting of the upper plane, so that practically all of the 
force accumulated in the stick phase is used to initiate 
the slip. 

Finally, the observed rate of force increase in the stick 
phase of the rotational regime (c) is much smaller than in 
the sliding regimes (a) and (b) for a similar reason. The 
rate of force increase in the stick phase is determined by 
an effective spring constant, K e ff, which is given by a 
combination of the elasticity of the spring n attached to 
the upper plane and the spring constant of the nanopar- 
ticle's contacts with both planes, K con t- Since this "con- 
tact" spring is attached to the spring of the upper plane 
in series, the combination rule is 1/ K e ff = 1/k+ l/ K eont, 
meaning that n e ff < K cont . In view of the large difference 
in the normal load, the contact in the sliding regimes (a) 
and (b) is much more rigid than in the rolling regime (c) . 
This implies a much smaller effective spring constant in 
the case (c) than in the cases (a) and (b), and a smaller 
rate of force increase. 

Fig. Q] shows that in the rolling state, the shape of 
the force curve [Fig. |3)Ja)] can be quite different from 
the saw-tooth-like ones shown in Fig. |3l Here, the ro- 
tation angle [Fig. HJb)] in the stick phase has the value 
<j> = mr/2 + 7t/4, meaning that the nanoparticle contacts 
the planes along its "incommensurate" facets. In the 
end of such a stick phase, the nanoparticle first quickly 
rotates by an angle of 7r/4 and enters another short- 
lived stick phase, where the contact is formed along the 
"commensurate" facets of the nanoparticle. Then, an- 
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FIG. 4: Temporal evolution of (a) the friction force / from 
Eq. ([3} and (b) the rotational angle <j> in the rolling state 
corresponding to the pulling velocity V = 4 and normal load 
In =30. 
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FIG. 5: Evolution of the friction force / from Eq. ([2j) during 
switching of the friction regime: (a) switching from HFS to 
LFS is performed by means of applying a negative pulse of 
the normal load at a constant velocity V — 2; (b) switching 
from LFS to HFS is achieved by applying a positive velocity 
pulse at a constant normal load /jv = 100. 



other rotation into an incommensurate state occurs, and 
the process repeats itself. As a result, the force evo- 
lution curve acquires an additional structure, with dif- 
ferent stick phases corresponding to different particle- 
surface contacts. We note that the force curve observed 
in the carbon nanotube rolling experiments— also pos- 
sesses a rather complex structure, presumably because 
of the variation of the contact properties of the nanotube 
"facets" . 



D. Control of friction 

The coexistence of the LFS and HFS states in Fig. [5] 
opens the possibility of switching between them by 
changing the control parameters, /jv and V. This pro- 
cess is illustrated in Fig. [5) Starting with the HFS-state 
within the coexistence region at fn — 100 and V = 2, 
we reduce the normal load to Jn — 30 keeping the ve- 
locity constant, Fig. EJa). This brings the system into 
that region of the state diagram from Fig. [2] where the 
LFS-state is the only stable one and induces a rotation 
of the nanoparticle by ir/4. Restoring the normal load 
to its initial value brings the nanoparticle back into the 



coexistence region, but now its contact with both planes 
goes along the incommensurate sides. 

In order to switch the nanoparticle back into the HFS 
state, one can increase the pulling velocity to a larger 
value, thus bringing the system out of the coexistence re- 
gion into the HFS part of the state diagram, see Fig.[5jb). 
This results in a sudden stretching of the spring and, cor- 
respondingly, in a large spike of the elastic force. After 
the nanoparticles has rotated into the HFS configuration, 
the velocity is reduced to the initial value in the coexis- 
tence region of the state diagram. Again, this velocity re- 
duction results in a sudden relaxation of the spring and 
in the large negative spike in the friction force. After 
that spike, the friction force stabilizes at a high value 
corresponding to the HFS-state of the system. 



IV. CONCLUSIONS 

To summarize, we have considered the rolling and slid- 
ing motion regimes of a cylindrical crystalline nanopar- 
ticle between two crystalline planes and found that the 
rolling state of motion is stable for sufficiently weak nor- 
mal loads and fast pulling. An intriguing feature of the 
state diagram from Fig. [5] is the existence of the stability 
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island of the LFS state, where the contact between the 
nanoparticle and the planes is incommensurate. This is 
in striking contrast with the finding of Refs. Hfla . where 
the HFS commensurable state was the only stable one. 
The reason for this difference is that, in the works Il4lfl5l . 
the flat nanoparticle was rotating around the axis perpen- 
dicular to the planes, while in our work, the rotation axis 
is parallel to the planes and perpendicular to the direc- 
tion of motion. We have shown that using nanoparticles 
as a lubricant, one can achieve either low or high friction 
at the same values of normal load and pulling velocity. 
This finding is at variance with the results from Refs. y- 
y, where adjustment of these parameters was essential for 
friction control, so that different friction regimes could be 
realized only for different values of V and /at. 

Even though the results reported in this paper have 
been obtained for the special case of equal lattice con- 
stants of the nanoparticle and the planes, we have found 
qualitatively similar behaviour also in the case of un- 
equal lattice constants. The fact that some facets of the 
nanoparticle are commensurate or incommensurate with 
the surfaces facilitates the discussion of the results, but 
is not central for the physical mechanism underlying the 
friction states of the nanoparticle. What is important 
is that different facets of the nanoparticle are character- 
ized by different interaction energies with the surfaces. If 
this condition is fulfilled, then several sliding states and 
the rolling state can be realized, even when all facets are 
incommensurate with the surfaces, or if they are com- 
mensurate, but differ in the number of contact atoms. 

In our analysis, several important effects have been 
neglected, such the nanoparticle 's possible asymmetry, 
inertia effects, and the effect of thermal noise. All these 
effects, when properly taken into account, may lead 
to qualitatively new friction regimes. In particular, 
nanoparticles of other shapes can exhibit different state 
diagrams and, correspondingly, can allow for different 



friction switching mechanisms. For instance, if the 
nanoparticle is asymmetric, it can stick to one plane 
and slide against the other; or it can slide with respect 
to both planes, but with the velocity very different from 
half the velocity of the upper plane. Next, the effects of 
thermal noise are typically negligible compared to load 
and interaction forces. However, in the systems where 
noise effects play a significant role, the nanoparticle is 
expected to spontaneously perform thermally induced 
transitions between the LFS and the HFS states within 
the coexistence region, so that higher normal loads 
would be required to stabilize them. In this case, one 
can expect the appearance of new regions in the state 
diagram, where friction is controlled by thermal noise. 
Finally, inertia introduces new characteristic time scales 
into the problem - the inverse resonance frequencies 
of the nanoparticle 's vertical, horizontal, and angular 
oscillations within the potential of the two surfaces. If 
the time of pulling by one lattice constant, a/V, becomes 
comparable to any of these time scales, new friction 
regimes associated with the nanoparticle's resonant 
motion can be expected. Exploring new friction regimes 
related to these and possibly other factors can be an 
exciting subject for future research, both theoretical and 
experimental. 
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